arXiv:hep-ph/0211257vl 16 Nov 2002 


MZ-TH/02-26 
November 2002 


Extraction of the strong coupling constant and 
strange quark mass from semileptonic r decaysf 1 ] 


A. A. Pivovarov 

Institut fur Physik der Johannes-Gutenberg-Universitat, 

Staudinger Weg 7, 55099 Mainz, Germany 
and 

Institute for Nuclear Research of the Russian Academy of Science, 

117312 Moscow, Russia 

Abstract 

In this lecture I present a pedagogical introduction to the low-energy phenomenology of light 
flavors. The renormalization scheme freedom in defining QCD parameters is discussed. It is 
shown in some details how one can extract an accurate numerical value for the strong coupling 
constant from the r-lepton decay rate into hadrons. As a related topic I discuss some peculiarities 
of definition of the quark mass in theories with confinement and describe the strange quark mass 
determination from data on r-lepton decays employing contour resummation which is a modern 
technique of the precision analysis in perturbative QCD. 

1 Definition of QCD parameters a s and m s 

Quantum ChromoDynamics (QCD) as a field theory for describing strong interactions is given 
by the Lagrangian 

£qcd = X 9/(*7/^ + ga'fG^ - m f )q f - jG^G^ 
f 

where is a non-Abelian gluonic field and G° w is the field strength tensor [|1|]. There are six 
quark flavors q = (u,d,s,c,b,t), three of which ( u,d,s ) are called light while the other three 
are heavy. There is a close analogy with QED - the Abelian gauge theory for describing the 
electromagnetic interaction of charged leptons. The QED Lagrangian for charged leptons reads 

£qed = X K ir r ^ 1 + ey- mi)l - jF^F^ 

i 4 

where F fW is the electromagnetic field strength tensor. In the standard model there are three 
charged leptons l = (e, fi , r): electron and muon are very light at the hadronic mass scale of order 
1 GeV while the r lepton is rather heavy with a mass M r = 1.777 GeV Q. 

The interaction is given by a vertex in the Lagrangian and normalized to a coupling constant 
g s at a tree level of perturbation theory (PT). A full theory (beyond the tree level) introduces 
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a dressed vertex that eventually determines a coupling constant after renormalization. For the 
renormalizable models of quantum field theory as QED and QCD the PT dressing is straightfor¬ 
ward and can in principle be done at any finite order of the expansion in the coupling constant. 
However, the dressing procedure is not unique because renormalization of loops can introduce 
some freedom in the choice of finite parameters through a particular definition of counterterms Q . 

Leptons can be detected as asymptotic states in the scattering processes that allows one 
to relate the coupling - for instance, the fine structure constant a - and the particle mass mi 
to observables very directly. Thus, the electron-photon scattering at low energies can be used 
to define the coupling constant (the fine structure constant, in fact) a = e 2 /47r through the 
Thompson cross section. This is a natural definition directly through a physical observable. 
Lepton mass is also directly related to experiment: it can be defined as a physical mass of the 
asymptotic state (a position of the pole of the lepton propagator). In QCD there is a phenomenon 
of confinement and no asymptotic states of quarks and/or gluons can be observed. Only the 
colorless hadrons appear as the asymptotic states. Therefore, definitions of the coupling constant 
and quark masses in QCD are less direct than in QED. At the theoretical level of the given 
Lagrangian they are very similar though: a vertex for the coupling and the propagator for a 
mass. To determine numerical values for the coupling and quark masses one should turn to 
experiment. As there is no possibility to measure these quantities directly one should specify 
the research area as the definitions of the parameters can be adopted to specific experiments 
(basically to a corresponding energy scale). In this lecture I will talk about r-lepton physics 
which is the area of low-energy hadron phenomenology: hadronic states have an energy E < M T . 
The particles which can be observed, for instance, in the process of e + e _ annihilation are p, w, 
ip, 7T7T. In r decays r —> v + hadrons one can in addition see 7r, ai(1260), K, A"*(892). 

The primary difficulty for extracting the QCD coupling from experimental data is that the 
qqg vertex cannot be “directly” measured. Indeed, a typical process (vertex) with hadrons at low 
energies is p —» tttt which is not directly expressed theoretically through the quark-gluon interac¬ 
tion vertex. Therefore, extracting a s (defined in terms of quark-gluon vertex in the Lagrangian) 
from the experimental quantity as, for instance, the p-meson decay width T(p —► 7T7r) is highly 
nontrivial. One can also see a propagation of the pion but not that of a quark. No mass shell for 
the quark (or gluon) is seen in the experiment. Thus, in QCD there is no preferable definition of 
parameters related to experiment. Then the only guidance for the choice of a particular definition 
of the QCD parameters in PT is the technical convenience (and also some general requirements 
such as gauge invariance) 0j|. 

Presently dimensional regularization is overwhelmingly used in many loop calculations. The 
renormalization procedure is usually a minimal one - subtraction of poles in £ = (D— 4)/2 where D 
is space-time dimension (an arbitrary complex number formally introduced for the regularization 
purposes). This procedure is quite abstract and remote from experimental quantities. Thus, 
formally, the MS-scheme coupling constant is defined by a ! f s (p) = Z^ s (p)af where af is a 
bare coupling constant. This definition of the renormalized coupling constant is not unique. In 
the momentum subtraction (MOM) scheme a^ [OM (/i) ~ rqqg(p? = p\ = p| = p 2 ). Since the 
definition is not unique the renormalization scheme freedom emerges. It is controlled by the 
renormalization group and can be conveniently parameterized by the coefficients (32,3,... of the 
/3-function [j|, |j] 

d 

p—ct s (p) = 2 (3(a s ), /3{a s ) = -/3 0 a 2 - (3ia 3 s - foot; - (3 3 a 5 s + 0(a 6 s ) 

and, therefore, any particularly defined coupling constant a s s ch generates an associated (3 sch (a s s ch )- 
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function. The same is true for the definition of the quark mass. In the MS-scheme one defines 
the renormalized mass as m MS (g) = Z^{p)m B with an associated 7 -function (pd / dp)m(p) = 
27 (a s )m(p). Functions /3(a s ) and 7 ( 07 ) are known up to four-loop approximation in PT j7[. The 
pole mass is difficult to define for light quarks since numerically even the strange quark is very 
light m s ~ Aqcd- 

In fact, a QCD coupling constant can also be defined beyond PT through physical observables. 
For light flavor phenomenology it can directly be defined through the cross section of e + e~ 
annihilation 


1 + a s (s) 


cr(e" 1 


hadrons) 


a(e + e~ —> pp) 

while for heavy quark physics the definition based on the heavy-quark static potential oty{q 2 ) 
Vqq(cf) can be useful ||]. 


2 Kinematics of semileptonic r decays 


The differential decay rate of the r lepton into an hadronic system H(s) with a total squared 
energy s 


da(r —> vH(s )) 
ds 


•171 + .* 


M?) 


Ml 


P(») 


is determined by the hadronic spectral density p(s) defined through the correlator of weak cur¬ 
rents. For the (ud) current jjf i x ) = ^ 7^(1 — 75 )d one finds 

p(s)ds 


* J( T J^ + (0))e igx dx = (, q^ - q 2 g flu )U had (q 2 


),n had (q 2 ) = 


s — q* 


(1) 


with p(s) ~ Im n had (s + iO), s = q 2 . The function II had ((5 2 ) with Q 2 = —q 2 is calculable in 
pQCD far from the physical cut as a series in the running coupling constant a s (Q 2 ). 

NonPT effects (power corrections) are included using OPE for the correlator (||) at small 
distances as x —> 0 in Euclidean domain (that corresponds to large Q 2 ) through phenomenological 
characteristics of the vacuum such as gluon and quark condensates ||. The lattice approximation 
for the evaluation of the correlator II liad ((3 2 ) beyond PT can also be used []l(]]. This is a basis for 
theoretical description of semileptonic r decays in QCD. 

Integrating the function II had (z) over a contour in the complex q 2 plane beyond the physical 
cut s > 0 one finds that for particular weight functions some integrals of the hadronic spectral 
density p(s) can be reliably computed in PT [11, |l2|. Indeed, due to Cauchy theorem one gets 


f n( 2 !)cfo = / p(s)ds . 
JC J cut 


Using the approximation II liad (z)| ze c ~ n PT (^)|^ e c which is well justified sufficiently far from 
the physical cut one obtains 


/ n had (£)cfe = f p(s)ds = (f II PT (z)dz 
Jc Jcut Jc 


i.e. the integral over the hadronic spectrum can be evaluated in pQCD. The total decay rate of 
the r lepton written in the form of an integral along the cut 


RtS =0 


r(r H S =pu) 

r(r —»• Ivv) 


/ cut 


1 Ml) 


1 + 


2s 

Ml 


p(s)ds 


is precisely the quantity that one can reliably compute in pQCD [ 13|. 
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3 PT analysis in QCD 


For technical reasons (no overall UV divergence) a derivative of the correlator II liad (Q 2 ) is often 
used for presenting results of PT evaluation 


D(Q 2 ) = 


d 


dQ‘ 


n 


had 


c q 2 ) = q 2 J 


p(s)ds 
(s + Q 2 ) 2 


The PT expansion for the D-function in terms of a s (Q ) = a s (Q)/TT reads D(Q 2 ) = 1 + a s (Q ) + 
k\a s (Q) 2 + /c 2 «s(Q) 3 + ksds^Q) 4 . In the MS-scheme 
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kl = ~2£ 9 C( 3) ; 


58057 779 . . . 75 >/ , 

= w-« 3) + t« 5 > 


with ((-function equal £(3) = 1.202... and ((5) = 1.037... The above numerical values for k 12 
summarize the results of three and four loop PT calculations [14]. Numerically, one finds D(Q 2 ) = 
1 + a s + 1.64a 2 + 6.37a(j + k^a*. Coefficient k% is known only partly [jD|. It is retained to obtain 
a feeling for the possible magnitude of the O(af) correction. The decay rate of the t lepton into 
nonstrange hadrons is written in the form 


d^rS =0 


T(t -> H s=0 v) 
r(r —»■ luv) 


= 3\V ud \ 2 (l + 5p + 5np ) • 


Here the first term is the parton model result, the second term 5p represents pQCD effects. 
NonPT effects are small, 5np ~ 0, in the factorization approximation for the four-quark vacuum 
condensates which is quite accurate |l3|, 16 j. 


The experimental result R^s=o = 3-492 ± 0.016 leads to 6 C p P = 0.203 ± 0.007 [jT?t| . In the 
MS-scheme the correction 5p is given by the series 


=a s + 5.2023a 2 + 26.366a d + (78.003 + fc 3 )a^ + 0(a p ) 


with a s taken at the scale p = M r . Usually one extracts a numerical value for a s {M T ) by 
treating the first three terms of the expression as an exact function - the cubic polynomial 
a s + 5.2023a 2 + 26.366af = 5 e p P . The solution reads tt af{M T ) = af{M r ) = 0.3404 ± 0.0073 e3; p. 
The error is due to the error of the input experimental value S C p P . It is difficult to estimate 
the theoretical uncertainty of the approximation for the (asymptotic) series given by the cubic 
polynomial (higher order terms). One criterion is the pattern of convergence of the series 

S e p p = 0.203 = 0.108 + 0.061 + 0.034 + ... 

The corrections provide a 100% change of the leading term. Another criterion is the order-by¬ 
order behavior of the extracted numerical value for the coupling constant. In consecutive orders 
of PT 

af{M T ) LO = 0.6377, af{M T ) NLO = 0.3882, af(M T ) NNLO = 0.3404 
that translates into a series for the coupling constant 

af{M T ) NNLO = 0.6377 - 0.2495 - 0.0478 


One can take a half(??) of the last term as an estimate of the theoretical uncertainty. It is 
only an indicative estimate. No rigorous justification can be given for such an assumption about 
the accuracy of the approximation without knowledge of the structure of the whole series. The 
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uncertainty obtained in such a way Aa| t (M r ) t / l = 0.0478/2 = 0.0239 0.0073 exp is much larger 
than that from experiment. This is a challenge for the theory: the accuracy of theoretical formulae 
cannot compete with experimental precision. Assuming this theoretical uncertainty one has 

af{M T ) NNLO = 0.3404 ± 0.0239^ ± 0.0073 exp 

Theory dominates the error. Still it is not the whole story. Now one can choose a different 
expansion parameter. The simplest way is to change the scale of the coupling along the RG 
trajectory M T — * 1 GeV. In terms of a s (l GeV) one finds a series a s (l)+2.615a s (l) 2 + 1.54a s (l ) 3 = 
5 C p P . The solution for the coupling constant is a s (l GeV) = 0.453. The convergence pattern for 
the correction 5 e p P is 

dp p = 0.203 = 0.144 + 0.054 + 0.005 
and for the numerical value of the coupling constant 

a s (l) = 0.453 = 0.638 - 0.177 - 0.008 . 


Should one conclude that now the accuracy is much better? What would be an invariant criterion 
for the precision of theoretical predictions obtained from PT, i.e. finite number of terms of 
asymptotic series? Thus, one sees that the renormalization scheme dependence can strongly 
obscure the heuristic evaluation of the accuracy of theoretical formulae in the absence of any 
information on the structure of the whole PT series. The final result for the standard reference 


value of the coupling constant normalized to the Z boson mass Mz reads [18] 


Oi s {Mz) T = 0.1184 ± 0.0007 ea: p ± 0.0006/jg mass ± 0.0010 th=truncation 

with a theoretical uncertainty that mainly comes from the truncation of the PT series. The world 
average value given by the Particle Data Group reads [2] 

a s {M z ) p V D G = 0.1172 ± 0.002 . 

To reduce a renormalization scheme dependence of the theoretical analysis one should use 
several observables simultaneously [19]. Such a possibility has recently emerged in study of r 
decays since experimental data on Cabibbo suppressed (S = 1) channel appeared |2(J. For 
strange hadrons H $=i (us part of the weak current: Cabibbo suppressed decays) the decay rate 
becomes 

D _ * 7Jg = \iy) _ OITA 12 / -| I £/ I £/ \ 

RtS =i — — 3|Vus| (1 + dp + d NP ) . 

The first term (“1”) is the parton model result, the second term 5' P gives pQCD effects. Small 
s-quark mass effects for Cabibbo suppressed part of the rate are taken in PT at the leading order 
in the ratio m 2 /M/ since rn s -C M r 


TTi 

5' P (a s ,m s ) = 5 P (a s ) + m (a s ) 

with 5p(a s ) being a correction in massless approximation for light quarks that is well justified for 


nonstrange decays (ud part) since u, cL quarks are very light indeed m u + = 14 MeV [21, 22]. 

The correlator of the weak charged strange current jp(x) = ^7^(1 — 75 )s with a finite s-quark 
mass is not transverse 


i I dxe lqx (Tj p (x)ji(0 )) = q p q u Tl q (q z ) + g pv U g (q 2 
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Retaining the first order term of expansion in the small ratio m 2 /q 2 one finds the m 2 correction 
to the invariant functions Ii q ^ g {q 2 ) 


mz 


U q (q 2 )=U(q 2 ) + 3^-U mq (q 


n g(q 2 ) = -q 2 U(q 2 ) + -m 2 U mg (q 2 


where n(q 2 ) is an invariant function for the mass zero case. The functions Ii qg {Q 2 ) are com¬ 
putable in QCD perturbation theory within operator product expansion for Q 2 —> oo. Thus, 
the experimental data on r lepton decays are theoretically described by three independent in¬ 
variant functions (form factors) which can be analyzed simultaneously that may help to reduce 
uncertainties introduced by the renormalization scheme freedom. 

In the actual analysis one can factor out the renormalization scheme freedom to large extent 
by introducing an effective scheme with definitions of effective quantities a, m 2 , m 2 through the 
relations [23] 


dQ 2 


-m s (M;)Q 


-Q 

2 —n 

dQ 2 ma 

m 2 (M 2 )U. 


n (Q 2 ) = 

(Q 2 ) = 
(Q 2 ) = 


mq 


1 + a(Q 2 ), 

m 2 (M 2 )C g (Q 2 ), 
m 2 (M 2 )C q (Q 2 ). 


(2) 


Here C q ^ g (Q 2 ) are coefficient functions of mass corrections. They are conveniently normalized 
by the requirement C q , g (M 2 ) = 1. In terms of the MS scheme quantities a s = a s (M 2 ) and 
m s = m s (M 2 ) the effective parameters in eq. (p[) read 


a(M 2 ) 

m 2 g (M 2 ) 

m 2 JM 2 ) 


QLc 


7r 


+ k\ 


+ k 2 


+ fa 


+ 0(a; 


5 a.. 


= m 2 s (M 2 )( 1 + ^ + k 

6 7r 


9 1 


+ kg'2 


+ 0(aj)) 


m l(M 2 )(l + - + kqi (—*) + k q2 +C , ( a s))- 


Numerical values for the coefficients fa, k q2 are unknown though their estimates within various 
intuitive approaches can be found in the literature [24], 

For confronting the experimental data obtained in r decays with theory of strong interactions 
one uses special integrals of the hadronic spectral density p(s) (called spectral moments) of the 
form 

rMl / „ \ k 


Mtl p = 


Pis) 


__s_Y 

M^J VM r 2 / M? 


which is suitable for experiment as the spectrum p(s) is measured. Another representation of the 
moments is suitable for the perturbation theory evaluation in QCD. Due to analytic properties 
of the functions n q :g (Q 2 ) the moments can be rewritten as contour integrals in the complex q 2 
plane: 

z \ k ( z \ l dz 

~Jd 2 ) KM 2 ) M 2 ' 

Technically one chooses a circular contour in the complex Q 2 -plane with Q 2 = M 2 e 1 ^, —ir < < ir 

that converts all invariant amplitudes Il q ^ g (Q 2 ) to the certain functions of the angle 0. The evo¬ 
lution of the functions H q ^ g (Q 2 ) along the contour in the complex plane is governed by the renor¬ 
malization group 125]. For the massless case corresponding to the analysis of data in Cabibbo 
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Figure 1: Running of the functions o( 4>) and C q (4>) on a circular contour in the complex plane 
calculated at LO, NLO and NNLO (left: a(</>); right: C g (^); real parts only) 


favored channel the only relevant quantity is the running “coupling constant” a(Q 2 ) —> a{4>) for 
Q 2 = M 2 e 1 ^ that serves as the expansion parameter for the function Yl{Q 2 ) along the contour. 
Note that there are no higher order corrections in the effective scheme by definition. The renor¬ 
malization group evolution is determined by the effective /5 function for the effective coupling 
constant [26, 27] 

= P{a(4 >)), a((j> = 0 ) = a(M 2 ), 

and the anomalous dimension for the running mass 

= 7 m[tt(^)]m s (^), m s ((p = 0) = m s . 

The initial values for a(</>) and m s (q 5) are extracted from fit to data. Thus, we have three quantities 
a(Q), C g (Q), C g (Q) associated with three invariant functions II (Q 2 ), H qt g(Q 2 ) describing the r 
system in the considered approximation. 

The RG equations for the set of quantities {a(Q),C q (Q),C g (Q)} are 

Q2 j^ a{Q2) = Q2 ^2 C ^ Q2) = 2 l9,q( a ) C g ,q(Q 2 ) • 

The RG functions /5(a) and 7 g, q (a) are given by the expressions [23] 

-4/5 (a) 


9a 2 


= 1 + 1.778a + 5.24a 2 + a d (-34 + 2 k 3 ); 


lg(yCl ^ = 1 + 4.03a + 17.45a 2 + a 3 (249.59 - k 3 ); 
a 

q — = 1 + 4.78a + 32.99a 2 + a 3 (-252 - k 3 + 3.4^). 
a 

The solution of the renormalization group equation for the effective coupling constant a(</>) con¬ 
verges well when the higher order corrections of the /5-function are included. The change from 
the next-to-leading order (NLO) solution to the next-to-next-to-leading order (NNLO) is small. 
The behavior of the coefficient function C g ((f)) (not shown) related to the contributions of spin 
one particles is rather similar to that of the coupling constant. However, the convergence pattern 
of the function C q (c/)) is much worse (see Fig. 1). It seems that the y^-function has already shown 
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Table 1: Coefficients of eq. (|3|) 


(M) 

aL U 
A kl 

4NLO 
A kl 

/tNNLO 
A kl 

n ki 

r nlo 

n kl 

rNNLU 

n kl 

(0,0) 

1.361 

1.445 

1.434 

0.523 

0.601 

0.625 

(1,0) 

1.568 

1.843 

1.976 

0.441 

0.552 

0.601 

(2,0) 

1.762 

2.282 

2.646 

0.390 

0.530 

0.607 


up an asymptotic growth in the next-to-next-to-leading order which will limit the precision of 
our results. Note that the fact that the function C q (4>) can behave wilder in higher order of 
perturbation theory is expected since this function is more infrared sensitive than the coupling 
constant a(4>) and the function C g ((f>). 

The numerical value for the s-quark mass m s is extracted from the difference between moments 
of Cabibbo-favored (ud- type) and and Cabibbo-suppressed (ws-type) decay rates 


8R™ 


R 


kl 

rS =0 


\v ud \* 


~nkl 

U rS= 1 i?kl 
IT/ 12 7 ^tS=0,1 

\ v us\ 



dR r s= 0,1 
ds 


The theoretical expression for the m? s corrections to the moments (k, l ) corresponding to the 
above experimental quantity is given by the contour integral in the complex q 2 plane 


3 im 2 s 

TV 


c 





II mg(z)\ dz 
M? ) M ?' 


In the theoretical expression for the difference <5/?/ we neglect terms of the order rnj /M 3 , set the 
u- and d-quark masses to zero, and retain only the most important term linear in m s . Within 
operator product expansion the coefficient of this term is given by the quark condensate. The 
final result for the difference reads |%|] 


m; 


5R k T ' = 3 Sew ( i^q^ki + u g B k {) - 47r 2 -^--^jT w 


(3) 


with S EW = 1.0194 |28] 


Here rrif i g = u q ^ g m\ 


with u) q = 1.73 ± 0.04, oj q = 1.42 ± 0.03. We use 


the relation between vacuum condensates of strange and nonstrange quarks 


{ss) = (0.8 ±0.2) (uu) 

and the numerical value (uu) = —(0.23 GeV) 3 [^]. In the leading order approximation of the 
QCD perturbation theory for the coefficient function of the quark condensate the quantities T k i 
multiplying the quark condensate are given by the expression 


T/d — 2 (Si,o(k + 2 ) - Si t i) . 


The numerical values for the first few coefficients T k i read 


2oo — 4, Tio — 6, T 20 — 8, Tqi — —2, Tn — —2. 


The numerical values for the coefficients A k i and B k i are given in Table 1. 
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Table 2: Experimental moments and extracted mass 


(M) 

(<ffif) ex p 

m s (Mf) MeV 

(0,0) 

(1,0) 

(2,0) 

0.394 ±0.137 

0.383 ± 0.078 
0.373 ± 0.054 

130 ± 5qq(= 6 ) 
111 ± <$$(=?) 
95±4(=?) 


Using experimental data one extracts m s as given in Table 2. Theoretical prediction for the 
moment (0,0) is the most reliable from PT point of view as <52o(=?) > ^io( = ?) > ^oo = 6 MeV. 
The final result reads 

m s (M T 2 ) = 130 ± 27 exp ± 3 (s - s> ± 6 th MeV. 

Normalization at 1 GeV gives 

m s ( 1 GeV) = 176 ± 37 exp ± 4^ s ) ± 9 th MeV . 

Only the moment (0, 0) is used for the m s determination as the most reliable one from the PT 
point of view. The higher order moments with the weight function (1 — s/M^) k for large k 
have an uncontrollable admixture of higher dimension condensates that makes them strongly 
nonperturbative and, therefore, unreliable for applications based on PT calculations |3(J. The 
contributions of higher dimension condensates are unknown and from general considerations the 
errors are expected to be much larger than 5 qq = 6 MeV. The value of the 

strange quark mass obtained by using the effective scheme approach as described in the present 
paper is in a reasonable agreement with other estimates |3l], |32|, |33|]. It is a bit larger than the 
recent lattice determination |]34j| . 

To conclude, the experimental information on r decays is a reliable source for the precision 
determination of the numerical values of important QCD parameters a s and m s . 
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